Compressed Sensing Matrices from Fourier Matrices 

Guangwu Xu *and Zhiqiang Xu ' 



Abstract 

The class of Fourier matrices is of special importance in compressed sensing (CS). 
This paper concerns deterministic construction of compressed sensing matrices from 
Fourier matrices. By using Katz' character sum estimation, we are able to design a 
deterministic procedure to select rows from a Fourier matrix to form a good com- 
pressed sensing matrix for sparse recovery. The sparsity bound in our construction is 
similar to that of binary CS matrices constructed by DeVore which greatly improves 
previous results for CS matrices from Fourier matrices. Our approach also provides 
more flexibilities in terms of the dimension of CS matrices. As a consequence, our 
construction yields an approximately mutually unbiased bases from Fourier matri- 
ces which is of particular interest to quantum information theory. This paper also 
contains a useful improvement to Katz' character sum estimation for quadratic ex- 
tensions, with an elementary and transparent proof. Some numerical examples are 
included. 

Keywords: l\ minimization, sparse recovery, mutual incoherence, compressed sensing 
matrices, deterministic construction, approximately mutually unbiased bases. 

1 Introduction 



In many practical situations the data of concern is sparse under suitable representations. 
Many problems turn to be computationally amenable under the sparsity assumption. As 
a notable example, it is now well understood that the l\ minimization method provides an 
effective way for reconstructing sparse signals in a variety of settings. 

The goal of compressed sensing (CS) is to recover a sparse vector (3 £ from the linear 
sampling y = $/3 and the sampling matrix (or compressed sensing matrix) $ £ C m . A 
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general model can be of the form y = <&/3 + z with z being a vector of errors (noise). 
In this case, one needs to approximate the sparse vector /3 from the linear sampling y, 
the compressed sensing matrix $ and some information of z. However we shall only be 
interested in the noiseless case (i.e., z = 0) in this paper to simplify discussion. 

In order to achieve less samples, one requires that the size of y be much smaller than the 
dimension of /3, namely m <C N. A naive approach for solving this problem is to consider 
£ minimization, i.e., 

(P ) (3 = argmin{||7|| subject to y — $7}. (1) 

However this is computationally infeasible. It is then natural to consider the method of 
i-y minimization which can be viewed as a convex relaxation of £q minimization. The i\ 
minimization method in this context is 

(Pi) (3 = argmin{||7||i subject to y — $7}. (2) 

7£<C iv 

This method has been successfully used as an effective way for reconstructing a sparse 
signal in many settings. 

Here, a central problem is to construct the compressed sensing matrices so that the 
solution to (Pi) is the same as that to (Pq). One of the most commonly used frameworks 
for the compressed sensing matrices $ is the restricted isometry property (RIP) which was 
introduced by Candes and Tao [8]. RIP has been used in the randomized construction of CS 
matrices (see for example [HEldD])- Another well-known framework in compressed sensing 
is the mutual incoherence property (MIP) of Donoho and Huo [H]. Several deterministic 
constructions of CS matrices are based on MIP, e.g., [121 I2U 137] . 

The main focus of this paper is on the deterministic construction of CS matrices. In [12] , 
DeVore presented a deterministic construction of CS matrices using the mutual incoherence: 
given a prime number p and an integer Kn^^aiimxiV (with m = p 2 ,N = p n ) binary 
matrix (i.e., each of its entry is either or 1) $ can be found in a deterministic manner 
such that whenever 



a signal /3 with the sparsity k (i.e., /3 has at most k nonzero components, we also call 
such a vector ^-sparse) can be produced by solving (PijE Recently, Li, Gao, Ge, and 
Zhang [24] suggested a deterministic construction of of binary CS matrices via algebraic 
curves over finite fields. As stated in [21], their construction is more flexible and slightly 
improves DeVore 's result when N is large (in fact the examples in [21] indicate that one 
sees improvement only when N > 0(m^")). 



1 In [12, it was stated that for RIP bound 6 < 1, it suffices that k < + 1. But for being able to 
recover fc-sparse signal, one needs to use © based on the discussion in [3]. 
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It is remarked that, in terms of construction via MIP, the bound is asymptotically 
optimal because of the Welch lower bound [M] (details will be given in section 2). This 
means that one needs to require k = 0(y/m). Using new estimates for sumsets in product 
sets and for exponential sums with the products of sets possessing special additive structure, 
Bourgain, Dilworth, Ford, Konyagin and Kutzarva [2] were able to overcome the natural 
barrier k = 0(y/m). However there is a restriction in the construction of [2], namely the 
ratio — is small. 

m 

Compared with binary matrices, the class of Fourier matrices is of great theoretical and 
practical relevance to compressed sensing. A notable technique of designing random partial 
Fourier matrices was suggested by Candes and Tao [9] , and was improved by Rudelson and 
Vershynin [29]. Let be the N x N Fourier matrix whose (k, j)-ih entry is given by 

For T = — Lj 7 ^), it was proved in [29] that if 

m = O (Hog 4 iV), (4) 

then a submatrix $ consists of m random rows of T satisfying RIP conditions (which 
ensures /c-sparse signal recovery) with high probability. This means that given N and m, 
the optimal k is k = O ( j^pr^ J ■ 

Deterministic construction of CS matrices from Fourier matrices has also received recent 
attention. In [35], by using the difference set, Xia, Zhou and Giannakis constructed a 
deterministic partial Fourier matrices with size m x N whose mutual incoherence constant 
meets the Welch lower bound (see Section 2). Such matrices are CS matrices for recovering 

a /c-sparse signal whenever k < \ ( \ ^jPr^ 1 + 1 J ■ However, the restriction on the matrix 
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size (m and N) is heavy. In fact, one can only construct such partial Fourier matrices for 
some very special parameters (m, N) with small — (see [35]). Recently, Haupt, Applebaum, 
and Nowak proposed a deterministic procedure to produce a partial Fourier matrix with 
size m x N for a large class of (m, N) pairs [19]. More specifically, the matrices in [19] are 
of size m x N with iV being a prime and m G [N 1 ^ , N] for some integer d > 2. This class 
of CS matrices can be used to recover fc-sparse signal with 

k = O (m^ 2 logd") . (5) 



It is noted that there is a significant gap between the bound ([5]) and the (asymptotically 
optimal) bound ()3]). The allowed range of sparsity for the case of a deterministic partial 
Fourier matrices is far smaller than that of the binary case. Constructing partial Fourier 
matrices deterministically that work for a larger sparsity range of signals is certainly of 
particular interest. 
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The aim of this paper is to construct a partial Fourier matrix which is a CS matrix, 
through a deterministic procedure. Based on a celebrated character sum estimation of Katz 
|20j . we are able to obtain a bound for the sparsity k in the case of partial Fourier matrices 
that is similar to ([3]). More precisely, we have shown that if q = p a is a prime power and 
n > 1, setting N = q n — 1 or N = = p p bZi {b I a is required for this case), then there 
is a deterministic process to select m = q rows from the Fourier matrix and build a 
(column normalized) matrix $, such that if 



Jm 1 . . 

then $ can be used to reconstruct A;-sparse signals via (-Pi). It is noted that ([6]) also greatly 
improves (jSJ). The result in this paper can also be used to recover the sparse trigonometric 
polynomial with a single variable j27J EE] • 

In this paper, we also improve Katz' estimation for quadratic extension fields with an 
elementary and transparent approach. Using this improvement, we are able to construct a 
CS matrix which is a partial Fourier matrix, and whose columns are a union of orthonormal 
bases. This is a useful construction for sparse representation of signals in a union of 
orthonormal bases which has been a topic of some studies (see, for example [T4" l [T5 | Wf \ ITS] ) . 
Moreover, this construction produces an approximately mutually unbiased bases which is 
of particular interest in quantum information theory. We also conduct some numerical 
experiments. The results show that the deterministic partial Fourier matrices has a better 
performance over the random partial Fourier matrices, provided that these two classes of 
matrices are of comparable sizes. 

This paper is organized as follows. The section below provides some necessary concepts 
and results to be used in our discussion. The main results are given in section [31 The 
discussion of computational issues and numerical results are contained in the last section. 



2 Background and Preparation 

2.1 MIP and Sparse Recovery 

Let $ be an m x N matrix with normalized column vectors $i, $2, • • • , Assume that 
each <3>j is of unit (Euclidean) length. The mutual incoherence constant (MIC) is defined 

as 

fi = max 

i>¥=3 
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Even though in many situations a small \x is desired, the following well-known result of 
Welch [33] indicates that /i is bounded below 



^ > 



N — m 
(N- l)m 



In [T4], Donoho and Huo gave a computationally verifiable condition on the parameter 
/j, that ensures the sparse signal recovery: let $ be a concatenation of two orthonormal 
matrices. Assume (3 is /c-sparse. If 

then the solution for (Pi) is exactly j3. 

This result of [2] was extended to a general matrix $ by Fuchs [TB], and, Gribonval 
and Nielsen [18], both in the noiseless case. For the noisy case of the bounded error, 
[TBI [51 [33] proved that ^-minimization gives a stable approximation of the signal (3 under 
some conditions that are stronger than (J7|). The open problem of whether condition (J7j) , 
namely \i < ni— r, i s sufficient for stable approximation of (3 in the noisy case was settled 
by Cai, Wang and Xu in [3]-actually the authors of [I] even proved that this condition is 
sharp, for both noisy and noiseless cases. Another remark is that if one considers only the 
noiseless case, the proof in [3] simplifies that of [HI [161 EH] 

Although the sparse recovery condition ([7]) is rather strong, it is advantageous in check- 
ing whether a matrix meets the condition. Such a checking procedure requires 0(N 2 ) steps 
which is computationally feasible. The MIC \x has been explicitly used in the design of 
compressed sensing matrices by several works, e.g., [3j [121 13 I2H [37] . 



2.2 RIP and Sparse Recovery 

We say that <3> satisfies the Restricted Isometry Property (RIP) of order k and constant 
5 k e [0, 1) if 

(l-^)Nl<lM<(l + (J fc )||a;|^ ( 8 ) 
holds for all fc-sparse vector x (see [S]). In fact, (JSJ) is equivalent to requiring that the 
Grammian matrices has all of its eigenvalues in [1 — 5k, 1 + 5k] for all T with \T\ < k, 

where $t is the sub matrix of $ whose columns are those with indexes in T. It has been 
shown that, under various conditions on RIP constant 5k, such as 5k < | (see [6]), one can 
recover the /c-sparse signal by solving (Pi). 



2.3 Katz' Characters Sums Estimation and Its Improvement 

Let ¥ q n be a finite field of order q n where q is a prime power. A multiplicative character x is 
a homomorphism from the multiplicative group (F* n , •) to S 1 where S 1 = {z £ C : ||z|| = 1}. 



5 



By a trivial character we mean the function that sends every element of F*„ to 1. 

Let g be a primitive root of ¥ q n (namely, g generates the multiplicative group F*„). For 
each element u G F*„, we define the discrete logarithm of u with respect to the base g to 
be the non-negative integer m ( m < q n — 1) such that g m = u, and we write log u 



m. 



Let N — q n — 1. Then a nontrivial multiplicative character of F*„ is of the form 



Xa{u) 



2iri log„ u 

1 — w 3 - 



where 1 < a < N — 1 . 

For a nontrivial multiplicative character, a celebrated theorem of Katz [20] concerns the 
magnitude of summation of the character values over a special coset of ¥ q (as an additive 
subgroup of ¥ q n). More precisely, Katz proved that for a nontrivial multiplicative character 
Xa of F*„ , and for an element a of ¥ qn with ¥ q ™ 



¥ q (a), there is an estimation 



a, 



te¥ a 



<(n~ l)y/q. 



(9) 



The method that Katz used in [20] to obtain the estimate ([9]) was a geometric one. In 
[25], Li presented an arithmetic proof of ([9]) using the Riemann hypothesis for the projective 
line over finite fields. It is interesting to note, in the case of quadratic extension of the 
fields, we can get a more precise estimation of Katz sum with an elementary and transparent 
proof. This result will be a useful tool in our construction of special CS matrices. Our 
improvement of Katz estimation is stated as 

Theorem 2.1 Suppose that Xa is a nontrivial multiplicative character of F* 2 . For an 
element a of¥ q 2 with ¥ q 2 = ¥ q (a), we have 



te¥ q 

5>a(* 



a, 



a) 



-i, 



if {q-l)\ a, 
if(q - 1) | a- 



(10) 



11^ 



*6F„ 



Proof. To this end, we first prove that 



t<=¥* 



q-l\a 
q — 1, q — 1 | a 



(12) 



Let g be a primitive root in ¥ q2 . For < m < q-2, the fact that (g m ^ +1 ))^~ l ) = g m ii 2 ~ l ) 
1 implies 

{^(g+i) . o < m < q _ 2} = F*, 
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or, equivalently 

(q+1) \\og g t < x=>- t G F*. 

Therefore, (O follows. 

Now assume that (q — 1) | a. We have 

= ^Xa(t)= ^ Xa(ti + t 2 a) 

teF* 2 ti,<2eF 9 ,(t 1 ,t 2 )^(o,o) 

= ^ *a(*l + *2Q!) = ^ Xa(0 + ^ Xa(* 2 ) ^ X«(*2 **1 

t£F* t 2 6F*tiSF 9 



(13) 



teF* 



t2£F* 



tieF, 



t£F* 



< 2 eF* 



t£F„ 



( 



1 + X Xa(t + a) ) X) Xa(*) = 1 + X + a) ) (g - 1), 

\ iSF 9 



*6F„ 



teF* 



and hence ^i 6 F„ Xo(* + at) = — 1. 

Next we consider the case (g — 1) j a and we shall prove (fTUj) . It is not difficult to check 

Q/ Q£ 

¥* for ti,t 2 G F 9 and ti 7^ t 2 . In fact, if = s G F 9 , then s / 1. 



that 



ti — a 



£2 - " 



= — — 7-i G F ? . This is absurd as a ^ F g . Therefore by ffTB"]) we see that 



s - 1 



(g + 1) f log ( 



ti — a 
to - a' 



provided £i,t 2 G F„ and t\ 7^ £2- Now 



i£F„ 



teF„ 



27rilog g (t— a) 
<j 2 -l 



2-7T1 log 9 (t 1 — a) 
Q 2 -l 



5>° 

l *2 6Fg 



-2iri log s (t 2 -a) 
q 2 -l 



q+ Yl e ° 



5 9 t 2 -a 



^-1 



g 2 -2 

2wifc i , „ 2nij(q+l) 



q 2 -2 

e 1 - 1 

fc=l,g+lffc 
</ 2 -2 



fc=l j=l 



0-2 

2,rifc , , „ 27rij(q + l) 



e 9 2 - 1 — > e" 1' 2 - 1 



k=0 



3=0 



which yields 



^2xa(t -a)\ = y/q. 



te¥„ 
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3 The Main Results and Construction Procedure 



In this section, we shall discuss construction of CS matrices by deterministically selecting 
set of rows from the N x N Fourier matrix. 

In the following discussion, we let F , Fx, ... , F N _i be the rows of the N x N Fourier 
matrix .F (JV) , i.e., 

j-(iV) = ^ 

Suppose that M is a subset of {0, 1, ... , N—l}, i.e., M = {m , mi, ... , m r } C {0, 1, . . . , N- 
1}. Then we can define the partial Fourier matrix associated with M as 



T 



(N) ._ 



M 



mi 



\Fm r / 



Let q be a prime power and let a G ¥ q n be such that 



F,n = ¥ q {a). 

Assume that g is a generator of the cyclic group F*„, we then have 

Theorem 3.1 Let q be a prime power, and n > 1 be a positive integer. Let N = q n — 1 
and 

M = {m = log g - a) : t G FJ. 
Then the q x N matrix $ = -j=J T m ) has MIC 



< 



n — l 



Proof. For < j, k < N — 1 and j ^ k, the inner product of the j'th and kth columns of 
the matrix $ is 



i 9-1 

X ^ — 27rjm r 'i 2ttA: % i 

Q r=0 

1 ^-^ (j / ,.) 27rlo M t - Q ) i 



q ^ 



r=0 



*6F„ 
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. . . 2-7T logg 

Since j 7^ /c, = e 77 defines a nontrivial multiplicative character of F*„. 



By the Katz estimation ([9]) we have 



fi = max = - 



5> 



feF„ 



< 



?7 



We would like to point out that this result was also mentioned in [2]. However, Theorem 
13.11 only produces a q x (g n — 1) matrix and hence the matrix size in Theorem 13.11 is quite 
restrictive. Note that in DeVore's construction [12], the matrix size is p 2 x p n where p is 
a prime. But Theorem 13. II cannot always be used to produce matrices of size Q(p 2 x p n ), 
e.g., when n is an odd number. 

We can have more flexibilities in choosing the size N of the Fourier matrix J 7 ^. This 
provides a larger set of parameters for constructing partial Fourier matrices for compressed 
sensing. More specifically, we have the following: 

Theorem 3.2 Let q = p a where p is a prime number. Suppose that b is a positive integer 
q n _ 1 



with b I a. Let N 
M = 



where n > 1 is a positive integer and let 

mod N) : t k G F„ k = 0, . . . , q - 1}. 

1 T~(N) 



p u — 1 

{m k = log (4 - a 



Then the MIC for the q x N matrix <J> = -^J^m satisfies 



n 



Proof. To this end, we first prove that #M = q, i.e., m k 7^ mj if k 7^ j. We assume that, 
for some k 7^ j, m k = mj holds. Then 



a 



log (tj — a) (mod TV). 



This implies that 



a 



tj-a 



fj 



dN 



(14) 



for some integer d. Obviously g dN 7^ 1 as t k 7^ tj. Noting that (g dN ) pb 1 = (g qn x ) = 1, 
we see that 



g dN e ¥ pb C ¥ q . 



Solving ( 1T4"|) for a, we get 



a 



g dN tj — t k 



.dN 



1 e F„ 
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which implies that F g (a) = ¥ q . This is impossible as n > 1. 

Now we consider the q x N matrix $. The inner product of the jth and kth columns 
of the matrix $ is 



. . 1 ^ V 27rjmi 2nkin\ 1 ^ \ 

= -}^ eNe N = - 2^ 



e (i-fc)( P "-i)^i 



me A/ 



rneAI 



(j-fc)(p"-l) 2 ' rl ° g g (t - n)i 



where iV = iV(jo fe — 1) = g n — 1. Since for < j, k < N and j ^ A;, 

x(«) = e w-*)^)^ 

defines a nontrivial multiplicative character of F*„. Again, by the Katz estimation fl9]), we 
arrive at 



H = max.\($ j ,Q k ) 



a. 



te¥ a 



< 



n — 1 



Remark 3.1 1. According to our discussion in the previous section, a k-sparse signal 
f3 in model (QJ) can be reconstructed via the l\-minimization (Pi) as long as fi < 



2/fc-l ' 



k < 



The matrix $ in Theorem \3.1\ and Theorem \3 '.2\ can be a CS matrix for recovering k 
sparse signals if 

Vq i 1 

2(n-l) 2' 

This bound is the same as the case of binary CS matrices obtained by DeVore JlW . 
Also, this bound improves the one for subsampling Fourier matrices in FTP]/ greatly. 

2. In terms of matrix dimension, our theorems provide more flexibilities. In fTM/ , the 
(binary) CS matrices can be of size p 2 x p r for any r > 2. We can get similar 
dimensions for partial Fourier CS matrices too. Taking q = p 2 in theorem \3.1\ we 
can get a p 2 x Q(p 2n ) matrix; letting a = 2 and b = 1 in theorem \ 3.2l we have a 
p 2 x 0(j> 2n_1 ) matrix. 



3. For the number N 



g n -l 



P an -i 

nf>- 



)U —, as pointed out in f^, the coherence of a q x N 
random matrix with i.i.d. Gaussian entries is about 2y / '^p- ~ 2y / ^ an ~^ logp . Hence, 



when logp > 



4(an— 6) 

smaller than that of random matrices. 



, MIC of the deterministic matrix $ given in Theorem \3.2\ is 
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Finally in this section, we restrict ourselves to the case of quadratic extension. We shall 
construct a partial Fourier matrix whose columns form a union of orthonomal bases, by 
using our improvement of Katz' estimation. The construction of such matrix is also raised 
in quantum information theory 

Let q be a prime power and let a be such that ¥ q 2 = ¥ g (a). As before, we let N = q 2 — 1 
and denote 

M = {m = log g -a) :te FJ. (15) 
It is easy to see that ^ M. We are able to state 



Theorem 3.3 Let 



1 



-.T 



(N) 



C (g+l)x( g 2 -l) 



y/q + T Mu{0} 

where M is defined in [To]) . For each j — 0, . . . , q — 2, set 

Tj = {j + k ■ (q - 1) : < k < q}. 

Then we have 

1. For any < j < q — 2, $jv is an orthogonal matrix. 

2. For k\ G Tj l and k 2 G Tj 2 with ji ^ j 2 , we have 



(16) 



q+ 1 



< 1(^,^)1 < 



q + l 



Proof. For any k±,k 2 G Tj with k\ ^ k 2 . Since (q — 1) | k± — k 2 , using fflTT) of theorem 12.11 
we get 



\ r=0 



27rfc^m r i 27rfc2 m r i 

e ^ e ^ 



^~ V r=0 



q + l 



(l + (-l)) = 0. 



This shows <I>7\ is an orthogonomal matrix for any < j < q — 2. 

We next consider the inner product of , <3>fc 2 with k\ and belonging to different 
sets Tj,j — 0, . . . , q — 2. We have 



1(^,^)1 



q + l 
1 



9-1 



q + l 



e — *r~ 

■r=0 



e ^ 



e ~ 



r=0 
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Since (q — 1) \ k\ — k 2 , using f TTU]) of theorem 12.11 we get 

q + 1 q+ 1 

■ 

Remark 3.2 T/ie concept of approximately mutually unbiased bases (AMUBs) (see 
[H| [Hj/J arises from quantum information theory. An mxN matrix <3> (with columns of unit 
length and with m \ N ) is a AMUBs if (i) the set of columns o/$ can be partitioned into — 

sets of m vectors, and each set is an orthonomal basis forC m ; (ii) \{^kii ®k 2 )\ — y^+ (l) 
holds if and $fc 2 are taken from different bases. By Taylor expansion, we have 



q + l V <? + 1 9 + 1 V 9+1 

5*0 Theorem \3. 3\ says that we have actually constructed an AMUBs. 



4 Computational Issues and Numerical Example 

In this section we will first discuss the issues of finding a primitive root g in W q n and of 
computing discrete logarithm with respect to g. We then explain some results from our 
numerical experiments. 



4.1 Computational Issues 

Finding a primitive root of a finite field is an interesting problem. In [32] , Shparlinski gave 
a deterministic algorithm which returns a primitive root g of F g n in time 0(q^ +£ ). It is 
noted that there is no polynomial time algorithm available for this problem even under the 
General Riemann Hypothesis unless q is prime and n < 2 (see Shoup [31]). The idea of 
Shparlinski's algorithm is to construct a subset M C ¥ q n such that 

1. M contains a primitive root; 

2. |M| = O(q^). 

To identify a primitive root in M one just uses the following fact: an element g G F*„ is a 

q n -l 

primitive root if and only if g~ ~ ^ £ for any prime factor I of q n — 1. 

Computing a discrete logarithm over a finite field is also of great importance because 
the hardness of this problem is a base of several well-known cryptosystems. The index cal- 
culus methods provide efficient ways of solving discrete logarithm problem over finite fields 
(see [26]). However, since in our setting, the size of the finite field involved is far smaller 
than that in the cryptographical setting, we can use the following algorithm of Shanks [30] 
to compute the discrete logarithm in F* n . Again, denote N = q n — 1 and let K = 
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Computing log 9 u for u G F*„ 

1 For all < j < K, compute gi and store (j, gi) in list L. 

2 Compute g~ K 

3 For i = to K - 1 

Compute u(g~ K y 

if u(g~ K y equals some g J in the list L, return log^w = iK + j. 
The computational costs are 0(y/N) both in time and space. 

4.2 Numerical Examples 

The purpose of the experiments is to provide some performance comparisons of the random 
sampling and the deterministic sampling of Fourier matrices. We use the method introduced 
in this paper to produce a deterministic partial Fourier matrix. To compare with the 
deterministic partial Fourier matrix, we generate the random partial Fourier matrix by 
choosing the same number of rows (as that in the deterministic case) from the Fourier 
matrix in a uniformly random manner. We reconstruct the Fourier coefficients by OMP 
(see [27J [36] ) and BP, respectively. For BP, we use the optimization tools of CVX [7J. 
We repeat the experiment 100 times for each sparsity k G {1, . . . , 20} and calculate the 
success rate in Example 1 and Example 2. In Example 3, we compare the maximum and 
minimum eigenvalue statistics of Gram matrices of varying sparsity k = \T\ for a 

deterministic sampling matrix and a random sampling matrix. It is interesting to note that 
all the numerical results show that the deterministic partial Fourier matrix has a better 
performance over the random partial Fourier matrix, regardless whether the sparsity if 
within the theoretical range (i.e. k < \ ^ + lj) or not. 

Example 1 Take q = 29 and N = q 2 - 1 = 840. Consider the field F 29 2 = ¥ 29 [x]/(x 2 + 2). 
Denote x = x + (x 2 + 2) and choose a = x, g = x + 2. We have 

Mi = {log 9 (t - a) : t e F 29 } 

= {465, 1, 494, 649, 47, 507, 758, 610, 835, 244, 67, 204, 588, 519, 

332, 808, 351, 672, 456, 683, 776, 275, 470, 562, 3, 103, 761, 466, 449}, 

By Theorem \3 . 1[ we can build a 29 x 840 partial Fourier matrix $ = -^J 7 ]^ . To compare, 
we generate the random partial Fourier matrix selecting 29 rows from the 840 x 840 Fourier 
matrix uniformly randomly. 

Figure 1 depicts the recovery results with using the recovery algorithm OMP and BP, 
respectively. 
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Figure 1: Numerical experiments for the comparison of the random and the deterministic 
29 x 840 partial Fourier matrix. The left graph corresponds to the success rates recovering 
by OMP, whereas the right one depicts the success rate of BP. 

Example 2 We take q = p = 19 and N = = p 2 + p + 1 = 381. Then we have 

M 2 = {\og g (t - a) mod iV : t G F 19 } 

= {192, 208, 162, 165, 160, 39, 154, 141, 245, 356, 304, 311, 223, 40, 321, 68, 118, 174, 249}, 

where F 19 3 = F 19 [a;]/(x 3 + x + 1), a = x, g = x 2 + 2x. By Theorem \ 3.£\ we can build a 
19 x 381 sub-Fourier matrix $ = -^J 7 ^. Similar to Example 1, we also compare the 
deterministic sampling and the random sampling and show the result in Figure 2. 




5 10 15 5 10 15 

Sparsily Sparsily 



Figure 2: Numerical experiments for the comparison of the random and the deterministic 
19 x 381 partial Fourier matrix. The left graph corresponds to the success rates recovered 
by OMP, whereas the right one depicts the success rate of BP. 

Example 3 The aim of the numerical experiment is to compare the maximum and min- 
imum eigenvalue statistics of Gram matrices of varying sparsity M = |T| for de- 
terministic sampling matrix and random sampling matrix. In fact, RIP of order k with 
constant Sk is equivalent to requiring that the Grammian matrices has all of its 
eigenvalues in [1 — 5k, 1 + Sk] for all T with |T| < k. For every value M , sets T are drawn 
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Figure 3: Eigenvalue statistics of Gram matrices for the deterministic sampling 

matrixes and random sampling matrixes. The left graph shows the result of 19 x 381 
partial Fourier matrix in Example 2, whereas the right one depicts the result of 29 x 840 
partial Fourier matrix in Example 1. 

uniformly random over all sets and the statistics are accumulated from 50,000 samples. 
Figure 3 shows the maximum and minimum eigenvalues of $J$ r for #T G {1, . . . , 20}. 
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